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ABSTRACT 

This paper describes the derivation of a series of equations 
capable of predicting the kinetic temperature of the aerospace environ- 
ment from 100 to 60,000 kilometers above the surface of the earth. The 
equations, which are programmed on a GE 225 computer, can be used to 
predict temperature-height profiles for any time through December 1992L. 



NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 



NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 


Technical Memorandum X-53104 


August 10, 1964 


AN AUTOMATED MODEL FOR PREDICTING THE KINETIC TEMPERATURE OF THE 
AEROSPACE ENVIRONMENT FROM 100 TO 60,000 KILOMETERS ABOVE 
THE SURFACE OF THE EARTH 

By 

Robert E. Smith 


AERO-ASTROPHYSICS OFFICE 
SPACE ENVIRONMENT GROUP 
AERO- AS TRODYNAMI CS LABORATORY 
RESEARCH AND DEVELOPMENT OPERATIONS 


DEFINITION OF SYMBOLS 


Symbol Definition 

ap geomagnetic activity index 


FT 

FTl 

FT2 

FZ 

FZl 


S F 

S Y 


SS 

Z 

z m 


correction factor for diurnal heating effect 

correction factor for variability of daily mean value 
of solar flux 

correction factor for variability of the monthly mean 
value of solar flux 

altitude variation of thermosphere lapse rate 

altitude variation of correction factor for the diurnal 
heating effect 

daily mean value of 10.7 centimeter solar radio noise 
flux (x 10~ 22 watts per meter 2 per cycle per second) 

monthly mean value of 10.7 centimeter solar radio noise 
flux (x 10“ 22 watts per meter 2 per cycle per second) 

yearly mean value of 10.7 centimeter solar radio noise 
flux (x 10- 22 watts per meter 2 per cycle per second) 

annual mean sunspot number 

Local Standard Time in hours and tenths 

altitude in kilometers 

altitude of thermopause in kilometers 
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TECHNICAL MEMORANDUM X-53104 


AN AUTOMATED MODEL FOR PREDICTING THE KINETIC TEMPERATURE OF THE 
AEROSPACE ENVIRONMENT FROM 100 TO 60,000 KILOMETERS ABOVE 
THE SURFACE OF THE EARTH 


SUMMARY 


A set of analytical expressions for predicting a time, solar flux, 
and geomagnetic activity dependent kinetic temperature-height profile of 
the upper atmosphere above the earth's surface is derived from the data 
presented by L. G. Jacchia in Smithsonian Institution Astrophysical 
Observatory Research in Space Science Special Report No. 150, dated 
22 April 1964 [1], Although the data sample is far from adequate and 
the concept of kinetic temperature in this altitude range is very nebu- 
lous , the data that can be provided by this automated predictive model 
is required to insure the orderly refinement of space vehicle design 
criteria. Values of kinetic temperature predicted by this model are 
compared in Figure 1 with the 1959 ARDC Standard and the 1962 U. S. 
Standard Atmosphere. 


I. INTRODUCTION 

In a previous paper [2], a model for predicting density-height pro- 
files in the aerospace environment from 200 to 60,000 kilometers above 
the earth's surface was developed through an analysis of current satel- 
lite observations. This is the second model of the series whose final 
model will be one which predicts the vertical distribution of the molecu- 
lar weight of the aerospace gas and its temporal and spatial variability. 
An approach similar to the one employed in Reference 2 was used in this 
development.. 


II. DISCUSSION 

The first step was the fitting of a curve by the least squares method 
to the data presented by Jacchia in Reference 1 combined with the kinetic 
temperature at 100 kilometers quoted in the 1963 Patrick Reference Atmosphere 
[3] assuming an almost isothermal structure above a thermopause whose 
height varies with solar activity [1], Equation (1), below, is a general 
equation for the temperature-height profile above 100 kilometers above 


the earth's surface at 0400 Local Standard Time (LST) when the 10.7 centi 
meter solar radio noise flux is 67 x 10” 22 watts per meter 2 per cycle per 
second and the geomagnetic activity index, ap, is zero. 


T 


(04,67,67,Z) 


(Z - 100)/ (0.002022168Z - 0.19844) + 210.00. 

( 1 ) 


Equation 2 will predict the height of the thermopause based on the 
monthly mean value of the 10.7 cm solar radio noise flux assuming a 
linear variation between Jacchia's heights of 220 kilometers at 
70 x 10" 22 w/m 2 -c-s and 400 kilometers at 220 x 10“ 22 w/m 2 -c-s [1], 


Z T = (6/5) S p + 135. 


( 2 ) 


Equation (3), below, which will predict a kinetic temperature 
envelope at any altitude, Z, at any time of the day, tL, for any solar 
radio noise flux, Sp and S F , and any geomagnetic activity condition, 
ap, incorporates all the corrective factors derived by Jacchia [1], 


T 


( t L»^D>^F , ^» a p^ 


T 


<04, 67 , 67, Z) 


+ (3.1 ± 1.9) a p 


1 + FT(FZl) 


+ FZ 


FT1 (S D - 67) + FT2 (S p - 67) . 


(3) 


The variability of the coefficient of ap is due to the latitude dependence 
of the geomagnetic activity effect. The lower value, 1.2, should be used 
in mid-latitude regions, 30° - 40°, while the higher value, 5.0, should 
be used in the auroral zones, 70° - 80°. If the ap values are unknown, 
then zero should be used in the equation and the resulting temperature 
will be the temperature predicted for a geomagnet ically quiet day [4], 


2 



In Equation (3), FT is a factor which will adjust the basic tempera- 
ture-height profile for the diurnal heating effect. The equation for 
this factor which varies from a minimum of 0.0 at 0400 LST to a maximum 
of 0.3 at 1400 LST is 



i0 + sin i0 


( 4 ) 


where 


a D = 0.1505 

a-L = 0.14347377 

a 2 = 0.8 x lO -2 

a 3 = -0.11666666 x 10” 2 

a 4 = 0.15 x lO" 2 

a 5 = -0.23071033 x 10” 2 


b x = 0.3674753 x 10" 1 
b 2 = -0.11835672 x 10” 1 
b 3 = -0.21666666 x 10” 2 
b 4 = -0.13333333 x 10” 2 
b 5 = 0.858005 x 10” 4 


a 6 = 0.000. 


In equation (3), FZ1 is the factor which makes the diurnal heating 
effect altitude dependent so that the diurnal variation increases from 
zero at 100 kilometers to 1.3 at and above 120 kilometers. 


FZ1 = (0.38 - 0.00275Z) (Z - 100.0) (5) 


for 100 S Z = 120 kilometers and FZl = FZl evaluated at 120 kilometers 
for Z > 120 kilometers. 
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In equation (3), FZ is a factor which makes the temperature lapse 
rate in the thermosphere altitude as well as solar activity dependent. 
The equation for this factor, which varies from 0.00 at 100 kilometers 
to approximately 0.85 at a 400 kilometer thermopause, is 


FZ = (Z - 100.0)/ (1.1096938Z - 92.828785) (6) 


for Z g Zp and FZ = FZ evaluated at Zj for Z > Zj. 

FT1, a factor which adjusts the predicted temperature for the daily 
variation, Sp, in the mean solar radio noise flux, varies from a minimum 
value of 1.9 at 0400 LST to a maximum value of 2.4 at 1400 LST. If the 
daily mean value of the 10.7 cm solar flux is unknown, then the value 
for Sp should be used. 

6 

FT1 = ^ c i cos + ^i s ^ n 
i=0 

where 


c Q = 2.15083333 

c-l = 0.23646679 
c 2 = 0.14166666 x 10- 1 
c 3 = 0.16666666 x 10" 2 
c 4 = -0.83333333 x 10" 3 
c 5 = -0.31334583 x 10“ 2 


d x = 0.58034166 x 10” 1 
d 2 = -0.24537375 x 10" 1 
d 3 = -0.16666666 x 1Q“ 2 
d 4 =* -0.16666666 x 10“ 2 
d 5 = 0.29916666 x 10“ 3 


c = 0.83333333 x 10~ 3 . 
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In Equation (3) FT2 is the factor which adjusts the predicted 
temperature for the variation in the monthly mean solar radio noise flux. 
The magnitude of this factor varies from a minimum of 4.5 at 0400 LST to 
a maximum of 6.0 at 1400 LSJT. 

6 

FT2 = ^ e i cos iQ + sin i0 (8) 

i=0 

where 

e Q = 5.26416666 

e x = 0.70135016 f x = 0. 16649887 

e 2 = 0.50833333 x 10" 1 f 2 = -0.11833333 

e 3 = -0.83333333 x 10" 2 f 3 = 0.00000000 

e 4 « -0.91666666 x 10 -2 f 4 = -0.14433705 x 10" 2 

e 5 = -0.30168333 x 10 -2 f 5 = 0.13501121 x 10 _1 

e 6 = 0.41666666" 2 . 

In equations (4), (7), and (8), 0 = (tL - 14.00) jt/12 where = Local 
Standard Time in hours and tenths on a twenty-four hour clock. 

The following set of equations was extracted from Reference 2. They 
will provide the solar activity predictions for use as input data. Equa-. 
tion (9) is used to convert either the known or predicted mean annual 
sunspot number, SS, to the mean annual value of 10.7 cm solar radio 
noise flux, Sy. Equations (10a) and (10b) then convert this mean annual 
value to the mean monthly value, Sj^, that would be predicted if the 
earth's orbit about the sun were circular. Equation (11) corrects this 
monthly mean value for the eccentricity of the earth's orbit according 
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to the data presented by Das Gupta [5]. These equations, when combined 
with a subjective prediction of the annual mean sunspot number from 
Figure 2, will provide all the data required to cotripute the predicted 
temperature-height profile envelopes. 

8 

Sy = Z 8 J SSJ (9) 

j=0 

where 

g o = 67.80036 g 5 = -0.261604186 x 10“ 8 

gi = 0.414498444 g 6 = -0.159036136 x 10 -1 ° 

g 2 = 0.427501471 x 10" 2 g 7 = 0.140461971 x 10 -12 


g3 = -0.469543503 x 10" 4 g 8 = -0.274461310 x 10" 15 


g 4 = 0.660862973 x 10" 6 


S 


M 


S Y + 



(S 


Y+l 


S^) when M § 6 


S 


M 




S Y ^) when H < 6 


S 


F 


M 


6 

^ tu cos i0 + sin i0 
i=0 


(10a) 


(10b) 


(U) 
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where 


h Q = 1.00 

h x = -0.28737186 x 10- 1 

h 2 = -0.10333333 x 10* 1 

h 3 = 0.150 x lO -2 

h 4 = -0.50 x 10- 3 

h 5 = 0.20601563 x 10" 2 

h 6 = 0.28333333 x 10” 2 

and 0 - 01-6) it/ 6 where M is the 
January = 1, February = 2, etc. 


k z = -0.13090493 x 10“ 
k 2 = -0.12701714 x 10* 
k 3 = -0.350 x 10' 2 
k 4 = -0.20207273 x 10“ 
k 5 = -0.24095063 x 10“ 

month of the year beginning with 


l 


l 


2 


2 


EDP Routine 

The following routine is programmed for use on a GE 225 computer. 
Input data for this program consists of the time of day, t^, the month 
of the year, M, the geomagnetic index, ap, the sunspot number, SS, the 
daily mean 10.7 cm solar radio noise flux, Sd, and the altitude, Z. The 
output consists of the time of day, TL, the geomagnetic index, AP, the 
sunspot number, SS, the month of the year, M, the altitude Z, the com- 
puted temperature, T, and the 10.7 cm solar radio noise flux, SF. 
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C MOD EL F 0 R cQdPUTl N 6 U P PER At HQ SPHERE TpHPERATURE 


fi B.E. S MITH JQ 6 NUMBER 55 0120 ._ 

COMMON SS1301, FHJ12J* 1 f 50 J , SFH f 7 1 , 8FK | 7 J > 6 f 9 J , 

IF 1 7 1 , M 7 ]• C ( 7 1 , D t 7 1 » A { 7 1 * i l 7l *SY 1301 

P16s, 52359878 

A*. 1*505 ....... 

A { 2 1 • . 14347377 

A 1 3 1 a « 8E-.2. „ 

A(4]a-.ll666666E-2 

A 1 5 1 a . 15E-2 — - 

Al6|a-,2307l033E-2 

_ AI7]*Q. - 

paQ, 

ilZl.*^i67.1Z53 E - 1 

Rl3!a-,ll835671E-l 

B M J « - . 2 1 6 6 6 6 6 6 E - 2 - 

R!5J»-,13333333E-2 

B 1 6 ] a * 858 Q 05E-4 ‘ 

B l 7 1 aO « 

XJLZ, 1518313! a.. 

C l 2 J ■ • 23646679 

Cl 31 * * 14 1.8 6 6 .6.61 - 1 — — .... 

C (4 ) a. 16666666E-2 

C l 5 ) a* » 833333331-3 , ... 

C16) a-. 313345831-2 

Cl 7 J ■ i .833333331-3 

D*0 , 

- 0 12 J a, 5 80341 6 6 E *1 

0131a-. 245373751-1 

D 1 4 J a*. 166666661-2. 1 - — 

D l 5 ) ■- . 1666666E-2 

IIi6.U_i.299l6A6.6E-3 

D l 7 ] aO • 

E*5,?6416666 _ ._ , 

El2)a. 70135016 

El 3 J a t 50 83.3.3.3 3.E-1. __ 

E14 Ja-.83333333B-2 

EL5 J a- .i9l.6jL4666.E- 2 

Fl6|a-.30l68333E-2 

E l 71 a * 416 6 66 66 E - 2 — 

F«0, 

F 12 1 a, 166.49887 __ 

Fl31a-. 11833333 

EXlXmJlj 1 

F 15! B-.14433705E-2 

_ .FI 6 J a. 13 5 0 112 IE* 1 _ — - _ 

F l 7 1 aO . 

- JGI67 .8 0 0 36 

Gl2!a. 419498444 

B. 1 . 3.1 m* 42ZS1L1.4 71E.-.2 , 

G14]b-,469543503E-4 

.615 J.a A .6 608 82973E-6 : 

GUI b*.261604186E-8 

6 171 a- .159 036136E-1Q - 

Gl81a,14046l971E-12 

G 1 9 1 a- « 2 7446 1311 - 15 , — 
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C MODEL FOft COMPUTING UPPER ATMOSPHERE TEMPERATURE 


>rn«i tu 

5 FH|?J«-.28737186E-1 
<?FH I 3 1 103333336*1 

SFH { 4 ) ■ .156*2 
sFHf 5l«-,5E-3 

SFHl6|«.2060l563E-2 
SFH 1 7 1 *. P83333333E*2 

SFK«0 • 

SFK 131 1309 0493E.1 

8FK131*-.12701714E-1 

SFKl4l*-.35E-2 

SFK 15 J •*.20207273E*2 
SFK 1 A 1 »*.24095Q63E*2 


A 

SFK ( 7 ) *0 . 
BCD, 10P 


5 

PRINT 5 
FORMAT l lMl 1 


A3 

60 TO 1 15*20 1 * I0P 

RCD.NM, NZ« TL, AP,$D,SS 111 . IFMl1 1 , I*1»NM| »(Z(I1*I *1* NZJ 



NS*1 

GO TO 30 


20 

RCD,NM,NS*TL» AP»SD,Z I 1 1 « iFM 1 1 1 # I«1#NM1 * |gS 1 1 1 # 1*1# NS I. 
NZ»1 

0 

SY CALCULATION * * ft ****** 

20 DO 35 1*1 j NS 

SYIM«fi . 

DO 33 K*g,9 


32 

SY 1 1 1 *SY 1 1 1 *G (K ) *SS( 1 1 ** I K»»l 1 


33 

CONTINUE 


40 

DO 130 IZ«1,NZ 

T1*I71 121*100. 1 / f 2. 0221AfiE*3*7 f 17 1-. 19844 1*91 0. 


41 - 

IFIZI IZ1-120. 141,41,42 

£RR«7 (12) 


42 

GO TO 43 
ERR«1 20 , 


43 

CONTINUE 

FZl«r.38-.00275*FRRl*[FPP-100. 1 



DO 145 IS*1*NS 
DO 140 IM«1,NM 



IFIFMI IM1*6, 180,50,50 


. 50 

JF t IS*NS170,60,60 


60 

SM«SYllSl 
GO TO 90 

C 

EQTN 

70 

9 A * * * * * * * * * * * * 

SM*SYl IS1 + IFMI IMl"6» l/12.*ISvl ! S*1 1 *SY [ IS) 1 


80 

GO TO 90 

I F C I s»l 160,60,85 


83 

.... 90 

SM*SVI ISl*t6.*FMl I H | )/12.*(SY(IS)-SY(iS*l) 1 
PHI«rFMHMJ*6. 1*PI6 



sf«sfh 



DO 95 1 «?, 7 



F I * 1 *1 
ANG*FI*PMI 


93 

SF*SFTSFHU |*COSFUNGl*SFKl 1 1*SINFIAN6I 

sf«sm*sf ; 



c MODEL FQ R COM p U T I N G UPPER ATMOSPHERE .TEMPERATURE- 

7T» 6.»S F/ 5.*135« ; 

TH» I TL”14 ♦ 1 *P 1 6/2 * 

FT»A ' __ 

FT2*E 

^ , FTi*r - m 

DO in5 I»2#7 

Filial - 

ang*f 1 *th 

SANGmSINF l ANG J . . 

CAnGbCOSFUNGJ 

FT2*FT2+E l I J*CANG+F1I1#SANS - - 

FTl*FTl+Cl IJ*CANG+Dm*SANG 

ins F T ■ F T ♦ A ! Li » CAN G ♦ J BjL 1 .l.t.SA MS 

iFIZf IZNZTJ110, 110*120 

110 FZMZl IZ1-100.]/ tU1096938*ZUZl-92. 8287851 

GO TO 130 

120 FZ«1ZT-100. 1/11.1096938 #xT - 9 2 • 8287851 * 

130 IFISD*1. 01131*131*132 

131 -Eliftasf : : _ 

GO TO 134 

132 ERRtSLi _ . - - - 

134 T«H.+FT*FZll*[Tl+l*2*APi+FZ*(FTt*[ERR-67.U 

IF 1 2* t SEhAZjJJ i -■ 

PH I NT 136«TL»AP»SS(1S)*FM(IM1*Z(!Z]« 

IXjlSF 

136 FORMAT t////2HTLEl7.9*2X2WAPEt7.9,2X2HSSEl7.9, 

12X2HM E17.9//2H Z El 7 .9 .2x2 HI E17.9*2X2HSFE17.V1 

140 CONTINUE 

145 CONTINUE — • 

150 CONTINUE 

G Q . T .Q..JL 

END 


IQ 


III. RESULTS 


Figure 3 depicts the diurnal variation of the temperature -height 
profile on a geomagnet ically quiet day in the mid-latitudes when the 
10.7 cm solar radio noise flux is 100 x 10“ 22 w/m 2 -c-s. 

Figure 4 depicts the diurnal variation of the temperature-height 
profile on a geomagnet ically quiet day in the mid-latitudes when the 
10.7 cm solar radio noise flux is 250 x 10” 22 w/m 2 -c-s . 

Figure 5 depicts the variation with solar activity of the tempera- 
ture height profile in the mid-latitudes at 0900 LST. 

Figure 6 depicts the variation with solar activity of the tempera- 
ture at mid-latitudes at 0900 LST at constant heights of 200 and 600 
kilometers above the earth's surface during the period from July 1964 
through December 1992. 
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AUitud* (x 10 2 fc») 

FIGURE 1. COMPARISON OF 1959 ARDC, 1962 U.S. STANDARD AND COMPARABLE 
TEMPERATURES PREDICTED BY THIS MODEL 
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FIGURE 2. PREDICTED MEAN ANNUAL SUNSPOT NUMBER FOR SOLAR CYCLES 20 AND 21 






T«apcrature (xlQ 2 °K) 


FIGURE 4. 


DIURNAL VARIATION IN TEMPERATURE -HEIGHT PROFILE AT 237 
UNITS OF 10.7 CM SOLAR RADIO NOISE FLUX 


x 10 -22 w/m 2 -c-s 
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FIGURE 6. VARIATION OF TEMPERATURE AT 200 AND 600 KILOMETERS DURING SOLAR CYCLES 20 AND 21 
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